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Abstract 

We calculate the supersymmetric 0{a s ) QCD corrections to the cross section e + e~ — > 
QiQj (hj = 1)2) within the Minimal Supersymmetric Standard Model. We pay particular 
attention to the case of the left-right squark mixing and to the renormalization of the 
mixing angle. The corrections due to gluino exchange turn out to be smaller than those 
due to gluon exchange, but they can be significant at higher energies even for a gluino 
mass of a few hundred GeV. 



1 Introduction 



Supersymmetry |l| requires the existence of two scalar partners (Jl and q~n (squarks) for every 
quark q. Quite generally, q~L and qn mix, the size of mixing being proportional to the mass 
of the quark ||. Therefore, the scalar partners of the top quark are expected to be strongly 
mixed so that one mass eigenstate ti can be rather light. Within the Minimal Supersymmetric 
Standard Model (MSSM) 6 £ and b R can considerably mix || |J for large tan/5 = v^/vx (where 
v\ and f 2 are the vacuum expectation values of the neutral components of the two Higgs dou- 
blets). Therefore, it is possible that ti or bi will be discovered in the energy range of the present 
colliders. 

The stop production in e + e~ annihilation, e + e~ — > t±ti, was first studied at tree level in [||]. The 
conventional QCD corrections to this process including virtual and real (soft and hard) gluon 
radiation were given in f| Q]. These corrections are quite large, for instance, they are about 
35% for m h = 80 GeV at y/s = 190 GeV §. The QCD corrections within the MSSM including 
virtual gluino and squark exchange were first calculated in |§, where the process e + e~ — > q~\q 2 
was also included. However, the renormalization of the squark mixing angle 9q as adopted there 
is not applicable in the whole range of 9q for the diagonal channels e + e~ — > qiq^ q2<l2- A proper 
renormalization of the squark mixing angle is necessary whenever stop and sbottom play a role. 
Moreover, no numerical analysis at all has been given so far for the unequal mass case or for 
squark mixing angle 9q^ 0. 

Here we want to present the complete formulae for the supersymmetric QCD corrections 
up to 0(a s ) within the MSSM including virtual and real gluon exchange, virtual gluinos and 
squarks for general (Jl^Qr mixing. In particular, we propose a suitable renormalization of the 
squark mixing angle. Furthermore, we give a detailed numerical analysis of these corrections. 
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2 Tree level formulae 



The squark mixing of qL and q R is expressed by 

?i = Ql cos 0g + q R sin 9g , g 2 = -q L sin 6^ + q R cos6»q , (1) 

where q±, q 2 are the mass eigenstates (with < mq 2 ) and 6q is the squark mixing angle. The 
production e + e~ — > g^-, (i,j = 1,2), proceeds via 7 and Z exchange in the s-channel (see 
Fig. la), s — (k + k) 2 , k and k being the four-momenta of the outgoing ^ and cjj. 
The cross section at tree-level is given by: 

o-°(e + e~ -> ft^.) = Cjj [eJ<Jij - T lZ e q a i:j 5 i: j + T zz a?-] , (2) 

with 
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TzZ 256 S 4 w i (a - m|) 2 + Y\m\ ' (5) 
Here \j = (1 — /i 2 — /i 2 ) 2 — 4/i 2 /i 2 , /i 2 ^ = m~. /s; e g is the charge of the squarks (e t = 2/3, e& = 
— 1/3) in units of e(= \J Aita) . a e and v e are the vector and axial vector couplings of the electron 
to the Z boson: v e — — 1 + 4s^ (with svk = siaOw), a e = — 1, and are the corresponding 
couplings Zqiq~j\ 

a u = A(P q L cos 2 0$ - s^ e<? ) , a 22 = 4(/| L sin 2 ^ - s 2 w e q ) , a 12 = a 21 = -llf sin 2# g ~ , (6) 
where is the third component of the weak isospin of the quark q. 

3 SUSY-QCD corrections 

The supersymmetric QCD corrected cross section in 0(a s ) corresponding to the Fig. 1 can be 
written as: 

a = a + Sa 9 + Sa 9 + Sa* . (7) 



Sa 9 gives the standard QCD gluonic correction (Figs, lb-f), 5a 9 is the correction due to the 
gluino exchange (Figs, lg, lh) and Sa 9 due to the squark exchange (Figs, li, lj). As will be 
seen later, within the renormalization prescription used Sa^ = 0. 
According to eq. (||D the correction can be written as: 



5a a = dj 2e q A(e q ) i i fS ij - T jZ (e q SijAa\j + A(e q )\fa ij ) + 2T zz a ij Aa 



a = g,g,q. 



3.1 Gluonic correction 



Sa 9 can be written as 



Aij has the following expression: 
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+ 2Li 2 (l - Aq) + Li 2 (A 2 ) - Li 2 (l - A 2 ) + Li 2 (A 2 ) - Li 2 (l - A 2 ) + 2 log 2 A - log A,, log A c 



with 



An 



(11) 



2/^i/^j ^Mi,j 

and Li 2 (x) = — Jq 1 log(l — dt. Eq. ( |i~0"|) has been calculated from the graphs in Figs, lb-f 
including soft and hard gluon radiation. For i = j the expression eq. ( |10D was already derived 
in H 0, |9|, |i~0"| . For « ^ j we also agree with the calculation of ref. apart from an obvious 
misprint in eq. (3.2) of 0. 
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3.2 The four— squark interaction 

Quite generally, we can decompose the corrections due to the four-squark interaction (Figs, li, 
lj) in the following way: 

AaSf = 8a^ + 8a^ + Sap, (12) 
A(e 9 )Sf = 5(e q )if > + 5(e q )\j^ . (13) 

Here the upper index v denotes the vertex correction (Fig. li), w the wave-function correction 
(Fig. lj), and 6 the renormalization of the mixing angle 6 q . The last term in eq. flT2] ) is 
necessary as the couplings a^- depend on the mixing angle (see eq. (|6|)). The vertex correction 
is proportional to the momentum of the exchanged photon or Z boson, and hence the related 
matrix element for e + e~ — > q^j is zero, that is 5ajj'^ = S(e q )^'^ = 0. There is, however, a 
squark wave-function correction 5a|*f due to the four-squark interaction. 
The self-energies due to the squark loops in Fig. lj are given by: 



la, / cos 2 26 q A° (m? ) + sin 2 26, A (m| ) \ sin 46 q ( A (m? ) - A (ml 



" ' 3 7t ^ |sin4^-(A°(m 2 - 2 ) - A°(m%)) sin 2 28 q A Q (m\) + cos 2 26~ q A*{m\ 2l 

' (14) 

Here A is the standard one-point function ||Tl |, A°(m 2 ) = —i7T~ 2 f d D q(q 2 — m 2 ) -1 in the 



convention of |T2|]. Notice that the self-energies are independent of p 2 , hence (m~J 



Sg^mD = Sjf. We therefore get for 5a^ ] : 

8a\j' q) = SZ^avj + bzf-a^ = -Re I l ' 1 a Vj + 2 J ] ' a if \ , ' J , (15) 



J ' .1 



m a ~ m a, m a ~ m a , 



with the squark wave-function renormalization constants 



r) Re {Eg) Reisf) 
^ff= / j , 2 l I ■ (16) 

Note that 8{e q )^ = because the contributions coming from the squark loop Fig. lj and the 
corresponding graph with the loop at the antisquark cancel each other. Hence A(e q )\f = 0. 
The correction term 5a in eq. ( |I2"D will be treated in section 3.4. 
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3.3 Correction due to gluino exchange 



Now we turn to the corrections due to gluino exchange (Fig. lg, lh). As in eq. (|TJ) and eq. (|l"3|) 
we can write: 



Aa 



(g) 



Salf* + 6a%' 9> + 5a\ 3 . 



S w >s) 



,0,9) 



A(e 9 )g> 



S(e q )[f } + 8{e q )% 



(17) 
(18) 



with 8a\j , 8(e q )^' 9> corresponding to the vertex correction (Fig. lg), and 5a \J' 9> , 8(e q )\J' 9> cor- 
responding to the wave-function correction (Fig. lh). Again 8afj is due to the renormalization 
of the mixing angle and will be calculated in section 3.4. 
The gluino vertex corrections are given by: 
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where 8^ is the identity matrix, v q = 2I q L — As^e^ a 
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The functions C^j are defined by 
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The arguments of all C-functions are (m~., s, m~ ,m~, m 2 , m q ). B°, C°, C 1 , and C 2 are the usual 



two- and three-point functions as given, for instance, in JT2 
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The squark wave-function renormalization due to the squark self-energy graphs with gluino 
exchange (Fig. lh) leads in the on-shell scheme to: 
1 



2 
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with the squark self-energy contributions £®(m 2 ) and their derivatives 
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3.4 The renormalization of the squark mixing angle 

We still have to discuss the renormalization of the squark mixing angle 5afj^\ eq. (12), and 
5a\j'**\ eq. (|iT|). From eq. (§) one gets for 5a (with a = q,g): 



5a\°' a) = 2a 19 5#> ) 



-5a 



5a 



(fl,a) 



5a 



(a 22 - an)56>, 



.(«) 



(27) 



Now Aa® of eq. (|12"D has to be free from ultra-violet divergencies and is therefore finite. 



Choosing Aa® = one gets from eqs. (0) and QlSf) 



5ag 5) = (a 22 - an W = -5Z^ J a 22 - 5Z^a n . 



.(«) 



(5). 



(28) 



Eq. means that the related counterterm 5a^ 2 '^ is nothing else than the negative sum of 
the graphs corresponding to Fig. lj containing the self-energies S® and £§? . We now do the 
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same for 8a[j by taking again the negative sum of those parts of the graphs Fig. lh, which 
contain the self-energies Sjf (m| 2 ) and E 2 f '(^J ■ 8a( 2 ^ is then also given by 

6af/ 9) = [a 22 - a u )W f ® = -8Z®a 22 - 8Z<§a n , (29) 

with 

r , Re|sg(m?)} ... ReiS^fm?)) 

2 ; , v 2 • (30) 

Inserting the results for 8Z[ 2 , 5^2?'* m ec L- 

(H) and 6Z$,6Z$ in eq. (|29|) , one gets 

S qM = f n4g g f^o (m a ) _ A o ( 2 j\ (31) 

' 6 7T m? -m| V v v , 



and 



S °f = lv i^(mf-ml) ( 5 °H.»^>n - *°K, m% m> 22 ) . (32) 

By eqs. ( p7|) we then obtain 5(4^ and <5a^'^' ) . With this result for 5aj^ there is no correction 
from the four-squark interaction graphs: Aa[f = 0, and due to A(e q )\f = also 5a^ = 0. 
By knowing 8a[j we can write the final result for the total gluino correction: 



Aa!ff = ^-^{^K^ (33) 
x ^(m? ,m?,mg)((-l) i+1 2a^cos2^- - ai/i/sin2^) + 5°(m|, , m~, mj)^ sin 2#J } , 



with i' ± i, and A(e 9 )g } = 8(e q ff + 8{e q )^ as in eq. (0), where <kjf \ S^(m|), 5(e ? )g' § ), 
and <5(e (? ). u ''^ are given by eqs. (||), (||), (0), and (H), respectively. Note that Aa\f and 



CU.C 

M A/zi ^s) 



A(e g )ij are ultra-violet finite. 



Inserting Aa-J and A(e 9 )^ in eq. (§) gives the total gluino correction 8a^ which does not 
factorize as the gluon correction in eq. @. 

Our renormalization condition for the mixing angle 6q is different from that of f9|. For 6g = 
and = rriq. our results agree with those of ||. However, in the scheme of [[J there appears 
a singularity at ^ = 45° going with ~ tan 26 q in the diagonal channels e + e~ — > <72<?2? so 
that this renormalization scheme is numerically not applicable in the region around 6q = 45°. 
In our scheme there is no such singularity as can be seen in eq. (^) so that it can be applied 
in the whole range of 6g. 
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4 Numerical results 



We now turn to the numerical analysis of the SUSY-QCD corrections. For the gluonic cor- 
rection Sa 9 we evaluate eqs. @, fllCf ) and (|TTD , and for the correction due to gluino exchange 



Sa 9 we take eqs. (|3^) and (|1^) together with eq. in the following we have always taken 
m t = 175 GeV. 

We have first calculated the corrections to the cross section a(e + e~ — > titi) at the LEP2 
energy ^Js = 190 GeV as a function of the stop mass for cos 6* t - = 0.7, taking m~ g = 200 GeV 
and mi 2 = 250 GeV. This is shown in Fig. 2. The gluon correction is rising from 17% of a tree 
for m~ tl = 45 GeV up to 35% for = 80 GeV. The gluino correction is only about 1% and 
quite independent of . 

In Fig. 3a we show the corrections to the cross section a(e + e~ — > tit\) as a function of the 
mixing angle cos# 4 ~, for a/s = 500 GeV, = 150 GeV, m t ~ 2 = 300 GeV, and m g = 300 GeV. 
According to eqs. (§) and (|T0|) 5a 9 has the same dependence on cos 9i as the tree-level cross 
section, whereas the gluino correction (see eqs. (|33|) ) introduces a different % dependence. This 
is due to the fact that the gluino couples equally to t* L t L and to t* R tR and therefore, in the case 
of mixing, the couplings gt*t get a dependence on the mixing angle % For the c. m. energy 
and masses chosen the gluino correction is now higher than at LEP2 energies and is about 
15% of the gluon correction. Fig. 3b exhibits the ^ dependence of 5a 9 /a tree and 5a 9 /a tree for 
cos#£ = 0.7 and keeping the masses as in Fig. 3a. Notice that 5cx 9 /a tree is decreasing with -^/i, 
whereas 5a 9 /a tree is becoming negative at = 700 GeV with the absolute value increasing 
with yfs. 

Moreover, we have examined the dependences of 5a 9 /a tree in e + e _ — > t\ti at fixed ^/s on the 
squark masses m^, mj 2 , and on their difference. These dependences are weak. 

Fig. 4a and b show the cos 9i dependence for a higher energy and mass scenario, y/s = 2 TeV, 
m h = 500 GeV, m h = 700 GeV, and m g = 600 GeV. Here 5a 9 is about -35% of 5a 9 for Q x 
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production. Fig. 4b exhibits the t\t 2 channel. 5a 9 /a tree is about 16%. 5a 9 is about —40% of 
Sa 9 at the peak values Notice that the gluino correction is not always positive. In Fig. 4c we 
also show for this scenario the corrections in the case of the ^2 production. Here the gluino 
part is even larger (—50% of the gluon contribution). 

We have also computed the SUSY-QCD corrections to a(e + e~ — > For ^/s = 2 TeV, 

mt = 500 GeV, = 700 GeV, and m g = 600 GeV 5a 9 is about —30% of Sa 9 almost inde- 
pendent of cos 6~ b . The gluon correction for bfij must be the same as that for t{tj production, 
provided the squark masses are the same. This can also be seen explicitely in eq. fllPp. 



It is particularly interesting to study the dependence on the gluino mass. This can be seen 
in Fig. 5, where we have \fs = 500 GeV, = 150 GeV, rri{ 2 = 300 GeV, and cos 6^- = 0.5 
in Fig. 5a, and = 2 TeV, m~ h = 500 GeV, m h = 700 GeV, and cos# t ~ = 0.5 in Fig. 5b. It 
is somewhat surprising that the gluino correction is first increasing as a function of the gluino 
mass and only very slowly decreasing. Of course the correction would be largest for a small 
gluino mass {m g < 50 GeV) but this is excluded by experiment {m g > 130 GeV). 
Here we also want to notice that threshold singularities appear when m- ti = m t + m g . For 
instance, such a singularity could appear in Fig. 5b for m g = 325 GeV. These singularities 
stem from the self-energy functions in the wave-function renormalization of the stops. This is 



a known effect and has been also observed in other cases [ 13 



In conclusion, our analysis of the SUSY-QCD corrections to scalar quark pair production 
in e + e~ annihilation has shown that the correction due to the gluino exchange are smaller than 
the conventional QCD corrections, but they are significant at energies envisaged for the next 
linear collider WM. In particular, they exhibit a dependence on the mixing angle 9g, which 
is different from the tree-level cross section corrected by the gluon exchange. Moreover, the 
correction due to gluino exchange decreases only very slowly with increasing gluino mass. 
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Figure Captions 

Fig. 1 Feynman diagrams for the lowest order SUSY-QCD corrections to e + e~ — > q^q^ Note 
that there are also the corresponding diagrams to c), d), e), h), and j) for the antisquark 

Fig. 2 SUSY-QCD corrections Sa 9 /a tree and 5a 9 /a tree as a function of m h for e+e" -> t Ji for 
v/£ = 190 GeV, cos# t ~ = 0.7, m t ~ 2 = 250 GeV, and m g = 200 GeV. 

Fig. 3 SUSY-QCD corrections 5a 9 /a tree and 5a 9 /a tree for e+e" -> i Jj 

function of cos#£ for yfs = 500 GeV, = 150 GeV, m ( - 2 = 300 GeV, and 
m g = 300 GeV, 

(b) as a function of y^i for cos6^ = 0.7, = 150 GeV, m,f 2 = 300 GeV, and m g = 
300 GeV. 

Fig. 4 SUSY-QCD corrections 5a 9 and <5a^ as a function of cos^ for 

v/i = 2 TeV, m fl = 500 GeV, m t ~ 2 = 700 GeV, and m g = 600 GeV. 

(a) for e + e~ — > titi 

(b) for e + e~ — > tit 2 

(c) for e + e~ — > t 2 ^2 

Fig. 5 Dependence of the SUSY-QCD corrections 5a 9 /a tree and 8a 9+9 /a tree on the gluino mass 
for e + e~ — > titi . 

(a) for ^ = 500 GeV, m h = 150 GeV, m h = 300 GeV, cos# t ~ = 0.5 

(b) for v/i = 2 TeV, m h = 500 GeV, m h = 700 GeV, cos 6~ t = 0.5 
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Fig. 1 
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